SISSA 1/98/FM 
physics/9801020 



Finite Mass Effect on Two Photon Processes in 
Hydrogenic Systems: Effective Scalar Photon 

Interaction 



C. K. Aii 1 ' 3 and C. S. Chu 2 ' 3 



1 Department of Physics and Astronomy, University of South Carolina, Columbia, SC 
29208, USA 

2 

International School for Advanced Studies (SISSA), Via Beirut 2, 34014 Trieste, Italy 
Department of Physics, The Chinese University of Hong Kong, Hong Kong, China 



email: au@sc.edu, cschu@tao.fm.sissa.it 



Abstract 



We consider a hydrogenic system with a nucleus of finite mass. The coupling of the 
radiation field to the center of mass motion gives rise to an effective scalar type coupling. 
This induced scalar photon interaction emerges as a correction in competition with the 
usual multipole interactions. This effect is particularly important in positronium where 
the electric quadrupole interaction is totally suppressed. We illustrate this effect with the 
two-photon decay of metastable hydrogenic systems. 



1 Introduction 



Recently, one of us |IJ pointed out, within the framework of nonrelativistic spinless QED, 
the emergence of an effective scalar photon interaction for the atomic (hydrogenic) internal 
degree of freedom due to finite mass corrections. In this simple model, the atom is a 
composite system with two particles of finite masses (mi and 777.2) and charges (Zi and 
Z 2 ) of opposite sign. In contrast to the relative motion which couples to the radiation 
field vectorially through minimal coupling, the coupling of the center of mass motion to 
the radiation field induces an effective scalar type coupling for the relative motion. Thus, 
when finite nuclear mass effect is taken into account, the induced scalar photon interaction 
leads to a modification of the multipole interactions of the atom with radiation. 

Consider a system consisting of two particles of masses mi, m 2 and charges Zie, Z 2 e. Let 
R and r denote the center of mass (cm.) and the relative coordinate vectors and P and 
p denote the momenta conjugate to R and r. To be definite, we refer to particle 1 as the 
nucleus and the particle 2 as the electron. The finite mass correction is characterized by 
the dimensionless quantity 

5 = m 2 /M, where M = mi + 7712 (1) 

is the total mass of the system. In general, < S < 1. The infinite nuclear mass 
limit corresponds to 5 = 0. It has been shown || that the spinless non-relativistic QED 
Hamiltonian describing this system is given by 

H = H + Hi nt , (2) 

where 

Ho = H 2part icle + H ra d, 

Hint = Hy + H s + H sg . (3) 

H is the Hamiltonian describing the two-particle system (H 2part i C i e ) and the radiation field 
(H rad ). H^t describes the interaction between the two-particle system and the radiation 
field. Hy describes the coupling of the radiation to the relative motion and is given by 
the usual vector photon type interaction. H sg IS db SCB. gull type interaction responsible 
for simultaneous two-photon processes. Hs describes the coupling of the radiation to the 
cm. and is given by an effective scalar photon interaction Hamiltonian. Explicitly, the 
forms of H v , H s , and H sg are [Q: 

Hy=Y, \[^ep ■ ee ^(^-D- _ ^L e ^ r) + h c (4) 

£t\ V U) 777 2 mi 
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where 




and e is the polarization vector and ci( e ,k) is the annhiliation operator for the photon. 

It has long been accepted that within nonrelativistic QED, a single photon transition 
between two S'-states (L = 0) is impossible M . This readily follows from parity consider- 
ation if only the vector interaction is assumed. In |]J , the effects of finite mass correction 
on one-photon transitions were studied. Taking into account the induced scalar interac- 
tion, it was shown that the metastable 2S state can decay by one-photon emission, which 
is otherwise forbidden if the nucleus is taken to be infinitely massive. Being a finite mass 
correction, this one photon decay rate is generally much smaller than the two photon 
decay rate. The Ml decay [[|, though allowed relativisticaly, is not considered here. 

For a self conjugate system with mi = m2 and Z\ = —Z2, it turns out that this finite 
mass correction effect vanishes in the one-photon decay channel and so the leading effect 
of finite mass correction will appear in two-photon processes. It is the aim of this report 
to provide a systematic study of finite mass effects on two photon processes. 

We now give an estimate of the role of the effective scalar photon interaction. For sim- 
plicity, we ignore the overall change in the state of the center of mass motion, bearing in 
mind that a nonvanishing effective scalar photon interaction requires a nonzero velocity 
for the center of mass motion. But in the calculations reported in sections 2 and 3 below, 
we do include the effects of the cm. recoil. As an example, we consider the two-photon 
decay of positronium from the metastable \2S > state to \1S > state. In this case, 
Z\ = —Z 2 = 1, mi — m 2 — m ei 5=1/2 and M = 2m e . On expanding the exponential 
factors exp(— ik-r/2) and exp(— ik-r/2) in both eqns. @ and (||), we readily observe that 
all even electric multipole interactions are totally suppressed in positronium M. We now 
focus on the dominant decay channels. For the El (electric dipole) channel arising from 
the vector interaction in the net result as compared to the massive hydrogenic case is 
that the electronic mass is replaced by the reduced mass in positronium (a difference by a 
factor of 2). Since the E2 channel is totally suppressed, the next contribution comes from 
the E3 channel giving rise to a correction to order a 4 (a being the fine structure constant) 
to the leading 2E1 decay amplitude. In contrast, on expanding the factor e ik r / 2 _ e * k r / 2 j n 
the effective scalar photon interaction in (|^), we see that the leading nonzero interaction 
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becomes 

H S ~H§> S - E J^.^-i e ^ (ik .r )a ^ + h ,.. (8) 

(e,k) ' 

We notice that as far as the atomic coordinate is concerned, the interaction in (|8|) behaves 
like an El (electric dipole) interaction whose coupling strength depends on the center of 
mass velocity. For comparsion, we also approximate the vector interaction Hy of ([|) by 
its leading El interaction: 

Hy~^ S -^ (9 ) 

Both interactions H^ and Hy^ couple to the \nP > state, but Hg has a nonzero 
coupling to the \2P > state whereas Hy does not because of the \2S > and \2P > 
degeneracy. Thus the \2P > state contributes to the two scalar photon decay of metastable 
postitronium, but not to the two vector photon decay nor to the one vector-one scalar 
photon decay. This follows from the trivial identity 

77?? 

< nP\p\2S >= ~Y^ n2 < nP\r\nS >, (10) 

where E n 2 is the energy difference between the \nP > and \2S > positronium states and 
we have explicitly used the fact that the reduced mass in positronium is m e /2. Thus we 
see that the contribution of Hjp to the two photon decay amplitude relative to that of 
Hy^ is of the order 

< H® > 1 kP P v 



< H$P > 4 m e E n2 m e c ' 

where v is the velocity of the positronium center of mass motion, and we have explicitly 
reinserted c. This simple analysis leads us to conclude that 

M ss M vs V 



M SV Myy C 



In this brief report, we consider a general hydrogenic system and study the effect of the 
scalar photon interaction H$ on the two-photon transition between two S 1 - states (nS to 
n'S). Thus, the general matrix element is of the form 

< n' S\H sg \nS > or 

< n S\HiGHj\nS >, where Hi(j) can be either Hy or H$, (13) 

and G is the hydrogenic Coluomb Green's function. For example, < n'S\H v GH v \nS > is 
the usual two vector photon transition matrix element between \nS > and \n'S >. This 



3 



can be calculated in closed form, following a method used by Gavrila and Costescu ||. 
More interesting are matrix elements of the type 



< n' S\H sGHy\nS >, the mixed vecor-scalar matrix element, 

< n' S\H s GH s \nS >, the two scalar photon transition matrix element. (14) 

The two scalar photon transtion matrix element for the special case where \nS > and 
\n'S > are both the \1S > ground state has in fact been evaluated analytically by one of 
us 0. It is interesting to point out that in the calculation reported in ||, the two scalar 
photon transition matrix element was only used as a calculational tool and no physical 
significance was attached to it. The generalization to arbitrary n and n' can be easily 
carried out and the result will be reported in section 2 below. The mixed scalar-vector two 
photon transition matrix element has never been studied before. It can also be evaluated 
in closed form and will be reported in section 2 below. In section 3, we apply these results 
to study the two-photon decay of the metastable 2S state. The angular distribution of 
the two-photon decay spectrum is reported. 



2 Exact two photon matrix elements 

We will see in the following that for a scattering process [] from (e, k) to (e', k'), it is useful 
to introduce the the following two photon matrix elements: 

the two vector photon matrix element, 

M vv (n, k 2 , k x ) = J J dridr 2 e' • pu n , s (r 2 ) e~^G{v 2 , n, ft)e ikiri e • pu^fa) 

dp 1 dp 2 e' • P2*Vs(P2 - k 2 ) G(p 2 , Pi, 0) e • Pi« n s(pi - ki), (15) 
the two scalar photon transition matrix element 

M 5S (fi,k 2 ,ki) = J J dp 1 dp 2 u n > s (p2 - k 2 ) G(p 2 ,pi,f2) u nS (pi - ki), (16) 
and the mixed scalar-vector two photon transition matrix element becomes 

M sv (£l, k 2 , ki) = / J dpidp 2 u n > s (p2 - k 2 ) G(p 2 , pi, Q) e • piw n5 (pi - ki), (17) 
M vs (Vt, k 2 , ki) = J J dp x dp 2 e' ■ p 2 w n / 5 (p 2 - k 2 ) G(p 2 , p x , ft) u nS {pi - ki). (18) 



1 Notice that for a two photon decay process, one only need to replace k by -k and uj by -u> in (|l5||l 
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Here G(j>2, Pi, is the Coulomb Green's function (Q — Hq) -1 in momentum representa- 
tion, ki = —5k or (1 — <5)k and k2 = 5k.' or —(1 — 5)k' originating from the finite mass 
effect. 

To simplify the discussion, we limt ourself to a hydrogenic atom where an electron is 
bound to a positive charge +Z\e\. We take advantage of the fact that the \nS > wave 
function can be generated from the \1S > wave function by a parametric differentiation: 

*-M = -^(^"(—^FW 1 ^' (19) 

where 

B n (x) = L l n (x)-2L 2 n (x), B x {x) = l, (20) 

and L v {x) are the associated Laguerre functions. In (|19"D, /3, n, X are independent variables 
and 

A = otZii. (21) 

Here /i is the reduced mass of the electron in the atom. We also use natural units where 
c = h = l. 

Next we used the fact || that the integral 

J(fl, k 2 , k x ) = J J d Pl dp 2 ((p 2 - k 2 ) 2 + r/ 2 )- 1 G(p 2 , Pl , fl) (( Pl - k x ) 2 + /3 2 )- 1 (22) 

can be evaluated analytically in closed form. The analytical form of I is given in (|3l| ) 
below. Then a comparison of fllTf ) and (122|), upon noticing the property of tranversality 



of the photon, clearly suggests that to generate the needed matrix elements, we need the 
simple relations: 

191 1 (23) 



2\2 ' 



2/3d/?((p + k) 2 + (3 2 ) ((p + k) 2 + /? 
Id 1 e ■ p 

~4 e ^((p + k) 2 + /3 2 ) = ((p + k) 2 + /? 2 ) 2 ' (24) 



Thus, the relations expressed in eqs. ( |23D and (|J) become the keys to generate a scalar 
photon vertex and a vector photon vertex respectively. 

Having laid down the rules to obtain the transition matrix elements, we proceed to simply 
stating the results. We define the following: 

X^t(2 f ,y/ 2 (E ns + Q) 1 / 2 , (25) 
E ns = -X 2 /(2fi)n 2 } (26) 

T EE X/X, (27) 
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P = X/n, T] = X/n', (28 
_ 2( 7 2 k 2 + (3 2 - X 2 )(k 2 2 + r] 2 - X 2 ) - 8X 2 kx • k 2 



(29) 



= ((X-/3) 2 + k 2 )((X-r ? ) 2 + k|) 
P ~ [(X + (3) 2 + k 2 ][(X + V ) 2 + k 2 } 1 ) 



and 

^ (1 — sp + pp 2 ) 



J(fi,k 2 ,kx) = 64A/iX(^— — ) J dpp~ T 
8A 5 



((X + /3) 2 + k 2 )((X + r ? ) 2 + k 2 ) 
^(^ka.kO, (31) 

where the integration contour begins at p = 1, runs along the real axis to a point closely 
on the right of p = 0, encirlces the origin in the counter-clockwise sense and runs back 
to p = 1. The integral J is expressible in terms of the Apelle type hypergeometric 
function. The desired matrix elements (p~5| ~ pT8|) are then expressible in terms of parametric 
differentiation of J. Specifically, 

M ™ = T^lkA, J(SlMM) ' (32) 

Mss = Wnhk J(nMMl (33) 
M - = 5^ £ '^ k - k ^ < 34 » 

M -4<^| J(n - k - k '»- (35 » 

This completes the analytical calculations of the two-photon transition matrix elements. 
We would like to point out that the parametric differentiation using eid/dki or e^d/dk^ is 
rather simple upon using the tranversality condition 

e . k = e' • k' = 0. (36) 



3 Two photon decay of the metastable 2S state 

In this section, we consider the two photon decay of the metastable 25' state of a hy- 
drogenic system with a electron bound to a positive charge Z\e\. We will present all the 
results consistently up to order A/M and (Za) 2 , compared to the leading order. 
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3.1 Kinematics 



Let the emitted photons be in the state (e, k) and (e',k'). Without loss of generality, we 
can take the atom to be initially at rest (P = 0). From the conservation of energy and 
momentum, we have 

p/2 

P' = -k - k' , (37) 

where A = ^fiZ 2 a 2 is the energy difference between the 25' and IS* level. This gives 

, ijj 2 + ijj' 2 + 2ujuj'cos6 . , - - . . 

u + J H — = A, where cos6 = k • k' (38) 

A convenient way to parametrize this equation is to introduce a dimensioness parameter 
a such that 

Acr(l - a) , A , o N1 

U = A[a + 2M (1 ~ C0S6) ~2M a + ° {a )] ' 

J = A[(l - a) + Aa ^ M a \ l - cos6) - A(i _ ff ) + o(a 8)]. (39) 

This parametrization satisfies ( |38| ) up to order A/M. The range of to is from to A — ^ 
and the range of a is from to 1. 



3.2 Transition amplitudes 

The transition operator is given by 



T = H 1 + H sg + H 1 - — #! + ••-, where H x = H v + H s . (40) 



1 

E — Hq 

It is clear that the matrix element for the terms in ellipsis in (|40|) are of order higher than 
(Za) 2 or A/M and we only need to compute < f\T\i > for the displayed terms in (flOf). 

The initial and final states are 

\i >= |P > ® \2S > ® |no photon >, 
|/ >= |P' > <g> \1S > ® I two photons in the modes (e, k), (e', k') > . (41) 

It is easy to see that 

</|ff 1 |i>=0, (42) 
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< f\H sg \i >= —{R(l -6) + Z 2 -^-R{5)}{e • ty^*** (43) 
in 2 1 — 

where 

R(5) =< lS\e- iS( ^ +k '> r \2S >, (44) 
which to order {Za) 2 can be written as 

R(5) = Ri(8) + R 2 (5)cos6, (45) 

where 

MS) = - °)(Za)H 2 ; (46) 

and 

< f\H 1 --^—H l \i >= e- l ( k+k ')- R (A(k,k') + B(k,k') + (k <-> k')), (47) 
£/ — N 

where A comes from the two vector photon interaction, B comes from the mixed vector- 
scalar photon interaction. They are: 

A(k,k') = -e 2 < 15|p ■ eY— e < (*- 1 > k '- r + — e i5k '- r )G(Q)(— e^' 1 ^ + — e i5kr )p ■ e\2S > 

m 2 m 1 m 2 mi 

£(k,k') = -e 2 ^-^ < 15|(Z 1 e l5k '- r + Z 2 e i(5 - 1)k ' r )G'(fi)(— e i(<5 - 1)kr + — e i<5k - r )p-e|25> 
M m 2 mi 

It is clear that these matrix elements can be reduced to the forms (|I5]) and (18). 
In particular, we have (for n = 2, n' = 1), 

Myy= 32 ^ X 3 [ - 0> + kf ) ^ " Ti 2 " ^ 3 - Ti *" *'> (4E 

/ 1 J2 2 — r 



F A (3-r;3,3;4-r;xi,X2) „ F A (4 - r; 3, 3; 5 - t;xi, x 2 ) u 

-Py* 3—7 ^ — T K e • 6 ) 

5^yV^ 2 2 yfl , 1 ^ ^(3-r;3,3;4-r;x 1 ,x 2 ) 
/I 1 /! ~ -A/^ + kJ ^— 

3 F A (4-r;4,4;5-r;a; 1 ,a; 2 ) 3 F A (5 - r; 4, 4; 6 - r; jci, £ 2 ) , w / , N 
~2^i 737 2^ 2 5^7 ](e ' 2)(e ' x) 

and 

64^A 4 X 3 „ . y 2P(X + P) ^ F A (2 - r;2,2;3- r;xi,x 2 ) 
Mvs= [ e ' k 2 ) [ 2{A + A)[L — ) 

, 4Qi/^F A (3 - r; 3, 3; 4 - r;xi,x 2 ) , , AQ 2 (3 ^ F A (A- r; 3, 3; 5 - r; Xi,x 2 ) 
/1 3-r /1 4-r 

3 gW^ F A (5 - r; 4, 4; 6 - r; xi, z 2 ) 3 ffWg F A (6 - r; 4, 4; 7 - r; x 1; x 2 ) 
2 A 5-r 2 /j 6-r J l49j 
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where Fa is the Apelle's type hypergeometric function, 

= A/2, A = aZp, (50) 

X 2 = 2p{\E 2s \+uj + ^) 1 (51) 

r = A/X, " (52) 

f 1 = (X + P) 2 + kl f 2 = (X + \) 2 + kj, (53) 

/ 1 = (X-/?) 2 + k 2 , f 2 = (X-\) 2 + k 2 2 , (54) 

g x = k\ + /? 2 - X 2 , ^ = k 2 + A 2 - X 2 , (55) 

2/i = 7 V[-4^ 2 X[(X + /3) 2 - k 2 ] + 16X 2 (X + /3)(ki ■ k 2 )], (56) 

Jl/2 

1/2 = tV[ 4X (( X " A ) 2 + W - X " ~ k ?)]' (57) 

Jl/2 

Qi = T^hWi/a + 3/3(X + A)/^ 2 + 3A(X + /3)f 29l - 5(X + (5){X + X) 9l g 2 

J 1/2 

-20X 2 (X + /3)(X + A)kx-k 2 ], (58) 
Q 2 = tV(- 4 ( X + - 4(X + \)PgJ 2 + 2\p gi g 2 + 6f~J 2 (X + f3)(X + A)] 

Jl/2 

-fi/ii4 (59) 

s = 7 ^-[2^ 2 + 8X 2 k 1 .k 2 ], P=£4 (60) 

Jl/2 J 1/2 

Wi = 2/12/2 + 2/22/i + sj/ij/J, W 2 = 2/22/2 - P2/i2/i (61) 

and Xi, x 2 are given by 

1 - sp + pp 2 = (1 - sip)(l - x 2 p). (62) 

Here y^, y' 2 are obtained from j/i, y 2 by interchanging A with /3 and ki with k 2 . 

The 2S to IS* amplitude Myy was first evaluated in [0] within the dipole approximation. 
Higher multipole contributions were calculated in ||. 



One can substitute these expression into ([40]) and obtain the exact two photon matrix 
elements. However, to illustrate the finite mass effect, it is enough to do an expansion 
and keep results up to first order in A/M. 

2 In ( jig] ) above, we have corrected some typos in ||. In ||, the total power of /1 in the first (second) 
factor of eqn. (7) should be 3 (4) instead of 2 (2); in eqn.(9), it should read x\ * x-i instead of x\ — £2; 
the L.H.S. of eqn. (13) should read t/2j the L.H.S. of eqn. (14) should read -y\. In ||, there was an error 
(a factor of 2) in one of the terms presented. The correct results are given in (^9h and (|72|) below. 
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It is convenient to write < 2S\T\1S > in the form 

<f\T\i >= ~^=M. (63) 

Using (|39D and perform a tedious but systematic and straight forward expansion in power 
of Za and A/M, one obtain 

M = (e- e'){g + (Za) 2 [f + (f + f"^)cos9]} + {Za)\h! + ^h")(e ■ k')(e' • k), (64) 

up to order (Za) 2 and first order in A/M. We have singled out the 5 dependence and 
g, f, f, f" and h! , h" are 5 independent. Here 9 is the correlation angle between the two 
photons in the initial rest frame of the decaying system. 

The physical origin of these terms are: g corresponds to the two electric dipole transition 
without retardation; / corresponds to correction to g to order (Za) 2 due to retardation 
effect; /', h! come from the two electric quadrupole transition. /' also contains a contribu- 
tion coming from the seagull term; /", h" include the recoil correction to the lowest order 
in A/M. 



3.3 Angular distribution of the decay spectrum 

On substituting mi by m e , the decay spectrum is given by 

dW = — ^ Y \M\ 2 uJu'dujdcos6. (65) 
2n mj 

Using (|39|) , one finds that 

uJdu = A 3 daa(l - <r)[l - -^-(1 + 2a) - -^(1 - a)cos6}. (66) 
Subsituting A = ^fi(Za) 2 = |(1 — S)m e (Za) 2 , one obtain 

(1U ^{l) 2 Z & m e a\l-5f- (67) 



dadcosO 2ir 8' 



Y \M\ 2 a(l - a) ■ [1 - — (1 + 2a)^-(Zaf - -(1 - a)cos6^-(Za 
' L 16 V ' M K 1 8 V 1 M K ' 
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where the sum is over the polarizations of the two emitted photons. To second order in 
Za and first order in A/M, we have 

y \M\ 2 = [g 2 + 2gf{Za) 2 }{\ + cos 2 6) + (Za) 2 [2g(f - h>) + 2g(f" - h")^]cos9 

+ {Za) 2 [2g{f + h>) + 2g(f" + h")^}cos*d. (68) 
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Hence, one can write 



^- = Z%d x + d 2 {Za) 2 + d^iZaf) • (1 + cos 2 9) 
dcosu M 

+{Zaf{d A + jjd 5 )cos6 + {Za) 2 (d 6 + j^d 7 )cos 3 6 } (69) 



with 



d\ = c J daa{\ — cr)g 2 , d 2 = c J daa(l — cr)2gf, 

d 3 = - c Jda3/16a(l - a){\ + 2a)g 2 , d 4 = cfda<r{l - a)2g(f - ti), 

d 5 = cjda[a(l - a)2g(f" - h") - 3/8(7(1 - a) 2 g% d 6 = cjdaa(l - *)2g(f + h% 
d 7 = cjda[a(l - a)2g(f" + h") - 3/8a(l - a) 2 g 2 ] 



and c = ^(§) 2 m e a: 8 ■ (1 - 5f 



(70) 



Integrating over a from to 1 by using the Simpson rule with 500 divisions and using 
2^(|) 2 m e a 8 = 52.3951 sec" 1 , we obtain (all numerical coefficients below are in units of 
sec -1 ) 

di = d i0 (l - 5)(1 - 5 + ZSf, for * = 1,3, 5, 7, 

d 2 = (1 - 5)(1 - 5 + Z5) 2 [((l - 5) 3 + Z 2 5 3 )a + (1-5 + Z5)((l - 5) 3 + Z5 3 )a% 
d 4 = d 40 {l -5)(l-8 + Z5) 2 ((l - 5) 3 + Z 2 5 3 ) 

d 6 = (l-8)(l-S+ Z5) 2 [((l - 5) 3 + Z 2 5 3 )b + ((1 - 5) 2 - Z5 2 )b% (71) 

with a = -4.3990 x 10 _1 , a' = -3.9877 x 10~ 2 b = -2.8271 x 10 _1 , b 1 = 7.2080 x 10~ 2 . 
The dependence on 5 in the coefficients d\ through d 7 comes from the recoil effects in the 
two-photons phase space (|66|) . In the infinite nuclear mass limit M = oo and 5=0, the 
coefficients g^'s become d i0 given by 

d w = 3.0860, d 20 = -4.7978 x 10~\ 

d 30 = -1.1573, d m = -2.8271 x 10 -1 , 

d 50 = 3.4852 x 10 _1 , d m = -2.1063 x 10" 1 , 
d 70 = -6.1464 x 10 -1 . (72) 



For positronium, Z = 1 and 5 = 1/2, the corresponding coefficients are given by 

d x = 1.5430, d 2 = -5.9972 x 10~ 2 , 
d 3 = -5.7865 x 10 -1 , d 4 = d G = -3.5339 x 10~ 2 , 
d 5 = 1.7426 x 10~ 2 , d 7 = -3.0732 x 10" 1 . (73) 

Notice that the electric quadrupole moment vanishes for the positronium and all the 
iZa) 2 corrections in ( |69|) arise from finite mass correction. It is not suprising that d A = d e 
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for the positronium if we remember that hi comes solely from the electric quadrupole 
transition and is zero in the present case. 
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